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[57] ABSTRACT 

Method and apparatus are provided for generating a simu- 
lation of a colored material. In particular, the method and 
apparatus generate a simulation of a colored target material 
utilizing a "selected" point in color space and a character 
reference pattern. The selected point in color space can itself 
be determined according to novel aspects of the invention 
(e.g., color candidate representation points), or other input- 
ted as a result of other processes (e.g., colorant formula 
performance analyses). In one mode of the invention, the 
selected point in color space can represent a color of a 
constituent element of the target material after the material 
has undergone one or more color-affecting treatments or 
processes. Advantageously, the present invention not only 
generates a simulation of a color material of a central color, 
but also generates simulations useful for representing color 
shade variations with respect to the central color. The 
method and apparatus are particularly suited for the simu- 
lation of colored fabric or textiles, such as denim, but are not 
limited thereto, having applicability to other colored mate- 
rials as well. 

45 Claims, 14 Drawing Sheets 
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METHOD AND APPARATUS FOR 
SIMULATING COLORED MATERIAL 

This is a continuation of application Ser. No. 08/227,642, 
filed Apr. 14, 1994, now U.S. Pat No. 5,493,518. 5 

BACKGROUND 

1. Field of Invention 

This invention pertains to method and apparatus for \q 
providing a predictive visible indication of ranges of shades 
for colored material, including ranges of shades within 
prescribed tolerances and ranges of shades statistically 
anticipated for a prescribed production specification. 

2. Related Art and Other Considerations is 
The concept of color space is useful in representing and 

modeling color phenomena. Color space is a three-dimen- 
sional space in which each point corresponds to a color, 
including both luminance and chrominance aspects. Hie 
tristimulus values R (red), G (green), and B (blue) form such 20 
a color space. As used herein, color space has three axes — 
L* (lightness); a * (red to green); and b* (blue to yellow). 
Tristimulus values are transformed (e.g., converted) into 
color space using, for example, the CIE 1976 L*a*b* 
equation. The CIE 1976 L*a*b* equation is explained in 25 
sundry prior publications including "CMC: Calculation of 
Small Color Differences For Acceptability", AATCC Tech- 
nical Manual/1992, pp. 322-324. 

In the production and/or use of colored materials, a color 
standard or specification is typically established for the 30 
material. It often occurs that the color standard can be 
achieved from a field of colorant formulations (e.g. a num- 
ber of variations of colorant mixtures). Modern computer 
color formulation systems (including spectrophotometric 
instruments coupled with computers) have been used to find 35 
concentrations of a set of colorants which, when mixed 
together, achieve such color standards. These modern com- 
puter formulation systems typically use the well-known 
Kubelka-Munk mixing laws. 

40 

Realistically, colorants themselves (e.g., colorant concen- 
trations) usually vary. Issues attending colorant variation are 
addressed in the thesis "Spectrophotometric Color Formu- 
lation Of Off Standard Colorants: Formula Selection With 
Respect to Match Sensitivity And The Optimization of 
Colorant Formulation With Respect To Colorant Availabil- 
ity", of Michael Keating in 1986. The Keating thesis pro- 
vides a model for optimizing an existing colorant formula- 
tion with respect to colorant availability and develops an 
index of "sensitivity", whereby a match of colorants can be 
chosen based on the relatively sensitivity of the match to 
changes in the concentrations of the colorants. The Keating 
thesis proposes an analytical method which determines 
parameters (e.g., shape, position, etc.) of an n-dimensional 
ellipse which describes sensitivity and optimized vectors of « 
concentrations. 

SheLyn Incorporated has developed a color control sys- 
tem which utilizes the Keating thesis. The SheLyn color 
control system predicts the performance (production pat- 
tern) of a given color formulation. As exemplified by FIG. so 
11, screen graphics and printed output of such a color control 
system provides a display of expected production variation 
of the formula against a CMC acceptability ellipsoid for that 
shade. 

The Keating thesis and SheLyn color control system deal 65 
only with colorant formulation, other factors such as stabil- 
ity and fading characteristics explicitly being beyond the 
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scope of the prior Keating research. Yet, in many industrial 
applications, subsequent treatment of colored materials 
affects not only lightness, but also produces a color shift. . 
Often a series of treatment steps or operations are required 
in route to a finished material. Generally a color standard is 
established with respect to the finished product 

One non-limiting example of such subsequent treatment 
occurs in textile production, particularly the manufacture of 
denim. As is well known, denim is typically produced using 
a dyed warp yarn extending in a first direction, and an 
uncolored fill yam which extends in a cross direction. Upon 
being woven with the warp yam and fill yarn, the unwashed 
denim is called "rigid" denim. The denim is then subjected 
to one or more washing operations at a laundry. The number 
and nature of the washing operations is dependent upon the 
desired commercial characteristics of the denim. The wash- 
ing operations can include (for example), in consecutive 
order, a rinse stage, a stonewash stage, an instant old stage, 
new age stage, and bleachout stage. 

A denim consumer (e.g., a blue jean manufacturer) gen- 
erally establishes a color standard for the consumer's prod- 
uct. For example, if a blue jean manufacturer produces a 
stonewash denim jean, the manufacturer will provide a color 
standard for the stonewash denim jean. 

For quality control purposes, the color standard is typi- 
cally in the form of a visual shade band or blanket. With 
respect to denim, the visual shade blanket is usually pro- 
duced in a manner such as the following: the laundry 
washes, to the pre-selected stage (e.g., stonewash, instant 
old, etc.), between one thousand and five thousand garments, 
cut from three or more rolls of sample denim. From these 
garments the laundry cuts and submits to the denim cus- 
tomer a range of swatches of differing coloration. The 
customer chooses a center swatch representing the preferred 
or ideal coloration, as well as other tolerance swatches 
which illustrate acceptable extremes of coloration. In some 
cases the swatches may then be mounted (e.g., in juxtapo- 
sition) on a substrate to form a blanket. 

After adoption of the color standard (e.g., production of 
the center and tolerance swatches), during denim production 
the color standard is used periodically for quality control at 
the laundry. For example, if, after reaching the desired 
washing stage at the laundry, a new batch of denim does not 
fit within the visual shade band represented by the blanket, 
the laundry must stop production and submit samples of the 
new batch to the denim consumer for approval. In view of 
such facts as variation between the shade range of produc- 
tion denim and the sample denim used for the shade blanket, 
as well as differences in concentration/formulation of dye 
for the warp yarn, the scenario of submission to the customer 
may repeat several times during the life of a product 
Unfortunately, each such stoppage occasions delay, as well 
as increased inventory and downtime at the laundry. 

The inefficiency of the current denim production/laundry 
quality control regime is exacerbated as customers demand 
not only higher quality denim, but denim of greater diversity 
(e.g., colorations diverse from conventional indigo, denim of 
differing character), all with shorter production lead times. 

Lacking in the prior art, among other things, is a system 
which provides the colored material (e.g., denim) consumer 
with a visible indication or prediction for the potential of 
long term shade variation before the colored material goes 
into production. 

Present techniques allow for color variation in order to 
achieve an image depicting physical characteristics or 
desired aesthetic modifications to a material design. For 
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example, U.S. Pat. No. 4,500,919 to Schreiber (incorporated 
by reference) discloses a system in which a color original is 
generated by taking into consideration aesthetically desired 
alterations. 

U.S. Pat. No. 4,954,976 to Noonan, entitled "Method For 5 
Simulating Dyed Fabric" (incorporated herein by reference), 
provides a simulated yiew (via color monitor or hardcopy) 
of dyed fabric. Noonan fails, among other things, to take into 
consideration variable shade ranges at differing stages of 
processing. Rather, Noonan uses color correction (e.g., 10 
shading) at edges of yarn shape to simulate the travel of yam 
over and beneath other yams. In particular, Noonan corrects 
colors of simulated yarn shape so that the average color is 
equal to the spectrophotometrically measured color of the 
actual yam shape. Thus, Noonan provides no suggestion for 15 
providing simulation of shade ranges of fabric, much less 
shade ranges at varying stages of processing. 



SUMMARY 



20 



Method and apparatus are provided for generating a 
simulation of a colored material. In particular, the method 
and apparatus generate a simulation of a colored target 
material utilizing a "selected" point in color space and a ^ 
character reference pattern. The selected point in color space 
can itself be determined according to novel aspects of the 
invention as discussed below (e.g., color candidate repre- 
sentation points), or other input as a result of other processes 
(e.g., colorant formula performance analyses). In one mode 3Q 
of the invention, the selected point in color space can 
represent a color of a constituent element of the target 
material after the material has undergone one or more 
color-affecting treatments or processes. Advantageously, the 
present invention not only generates a simulation of a color 35 
material of a central color, but also generates simulations 
useful for representing color shade variations with respect to 
the central color. The method and apparatus are particularly 
suited for the simulation of colored fabric or textiles, such as 
denim, but are not limited thereto, having applicability to 4Q 
other colored materials as well. 

In one mode of the invention, a calculation processor 
determines the "selected" point in color space corresponding 
to the desired color of the colored target material. In 
response to user or stored input in the form of shade 45 
variation tolerance (e.g., DE^c), an image processor deter- 
mines a three dimensional volume (e.g., an ellipse) of 
acceptable color variation about the selected point, as well 
as a plurality of shade variation representation points. The 
shade variation points exist in color space on a surface of the 50 
three dimensional volume. The shade variation representa- 
tion points together with the selected point form a set of 
color candidate representation points. The image processor 
drives an output device to generate a shade simulation 
output for each color candidate representation point, thereby 55 
forming a set of simulation outputs depicting a range of 
predetermined shade variations. The output device can be a 
printer or a display screen. 

As indicated above, in generating the shade simulation 
outputs), the image processor takes into consideration not 60 
only the selected point in color space (e.g., the desired 
color), but also the character of the type of material being 
simulated. As used herein, the term "character" refers to one 
or more physical properties of constituent element(s) of the 
target material, including sizing and positional relationships 65 
of the constituent elements. For denim, for example, char- 
acter refers to one or more of the following: size of warp 



yam, size of fill yarn, tightness of weave, pattern [e.g., twill, 
right hand or left hand]), lb obtain an indication of charac- 
ter, a character reference pattern is obtained either by 
analysis of a character reference material by a scanner, or by 
selection from a character reference pattern library stored in 
memory accessible to the image processor. The character 
reference pattern is stored in memory as an array, each of a 
plurality of elements corresponding to a sub-area of the 
simulation output (e.g., a pixel). Each element of the array 
is assigned a corresponding initial lightness value. Utiliza- 
tion of the character reference pattern, in conjunction with a 
particular one of the representation points in color space, can 
occur either in a linear or non-linear mode. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing and other objects, features, and advantages 
of the invention will be apparent from the following more 
particular description of preferred embodiments as illus- 
trated in the accompanying drawings in which reference 
characters refer to the same parts throughout the various 
views. The drawings are not necessarily to scale, emphasis 
instead being placed upon illustrating the principles of the 
invention. 

FIG. 1 is a schematic view of a shade simulation system 
according to an embodiment of the invention. 

FIG. 2 is a schematic view depicting input of swatches 
into a spectral analyzer of the embodiment of FIG. 1. 

FIG. 3 is a flowchart showing general steps executed in 
connection with a material simulation process utilizing at 
least one color candidate representation point. 

FIG. 3A is a flowchart showing general steps executed in 
connection with a material simulation process utilizing at 
least one formula performance representation point 

FIG. 4 is a graph showing spectral analysis information 
provided by a spectral analysis subsystem included in the 
embodiment of FIG. 1. 

FIG. 5A and 5B are graphs showing spectral analysis 
information for various washing stages of a swatch of fabric. 

FIG! 6 is a graph illustrating representation points in 
L*,a*,b* color space. 

FIG. 7 is a table showing numerical values of spectral data 
for one type of swatch. 

FIG. 8 is a flowchart showing steps included in the 
execution of swatch simulation generation. 

FIG. 9 is a diagrammatic view depicting re-justification of 
lightness coordinates of reduction shades about a lightness 
coordinate of a representation point according to a linear 
mode. 

FIG. 10A-FIG. 10E are diagrammatic views showing 
values assigned during differing steps of processing to 
elements of an array used to prepare a simulation of a 
colored target material. 

FIG. 11 is a diagrammatic view of output of a prior art 
color control system which provides a display of expected 
production variation of a formula against CMC acceptability 
ellipsoids. 

FIG. 12 is a flowchart showing steps included in a 
non-linear mode of execution of certain steps of swatch 
simulation generation. 

FIG. 13 is a graph illustrating an ellipse in only two 
dimensions (c, h) in a L-c-h coordinate system 

DETAILED DESCRIPTION OF THE DRAWINGS 

FIG. 1 shows a material simulation system which includes 
a spectral analyzer subsystem 20; a calculation processing 
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subsystem 22; a scanner subsystem 24; and an image 
processor subsystem 26. In some modes, the material simu- 
lation system of FIG. 1 operates in conjunction with a 
material treatment processor 28. 

As subsequentially explained in more detail with refer- 5 
ence to FIG. 2, special analyzer subsystem 20 (also known 
as analyzer 20) receives one or more colored materials for 
analysis and separately, for received each material, conducts 
a spectra analysis. Color material is, in the illustrated 
embodiment, typically received for analysis in the form of a 10 
swatch. As a result of the spectral analysis, for each received 
material, analyzer 20 generates data sufficient to generate the 
graph shown as FIG. 4. In particular, analyzer 20 provides 
a percent reflectance value (%R) for wavelengths in a range 
approximating the visible portion of the electromagnetic 
spectrum (e.g., from about 400 nm to about 720 nm). In 15 
addition, for each analyzed material, analyzer 20 provides 
values for L*. a*, b*. and C*, as well as an angle value 9, 
which are collectively referred to as spectral data. As well 
understood by those skilled in the an, L*, a* f b* are 
CEELAB standard nomenclature as explained above. 20 
Regarding other nomenclature, C* refers to Chroma and 
angle 6 refers to hue angle. In the illustrated embodiment, 
the spectral analyzer subsystem 20 is provided by Hunterlab 
and bears model number COS- 1600. 

25 

Arrow 30 in FIG. 1 represents input of spectral data from 
analyzer 20 to calculation processing subsystem 22. In one 
embodiment, input of spectral data from analyzer 20 to 
calculation processing subsystem 22 is by manual entry. In 
an automated embodiment, for each analyzed material, each 3Q 
of the values included in the spectral data are transmitted 
over suitable electrical cable(s) to calculation processing 
subsystem 22 in accordance with a predetermined protocol. 
For example, in response to a command signal from sub- 
system 22 carried on line 31, analyzer 20 conducts the « 
analysis of a swatch of analyzed material received therein 
and transmits spectral data for the swatch to subsystem 22 as 
depicted by arrow 30. 

Calculation processing subsystem 22 has connected 
thereto (via conventional interfaces) an input device 32, an 40 
output device 34, and memory 36. In the illustrated embodi- 
ment, input device 32 is a keyboard; output device 34 is a 
printer. As explained hereinafter, calculation processing sub- 
system 22 is a computer configured to control the sequence 
of operations of analyzer 20, image processor subsystem 26, 45 
and scanner 24, as well as to control transmission of data 
between itself, analyzer 20, and image processor subsystem 
26. Bus 38 carries signals and data between calculation 
processing subsystem 22 and scanner subsystem 24 (also 
abbreviated herein as "scanner'*). Line 39 is employed in the 50 
illustrated embodiment to control scanner 24. Calculation 
processing subsystem 22 also makes certain calculations as 
hereinafter described. 

As explained in more detail below, scanner subsystem 24 
receives a character reference swatch and, upon either 55 
manual or automated command (from subsystem 22 on line 
39), generates a two-dimensional color image or color map 
of the scanned character reference swatch. The color map is 
used as a color reference partem. That is, scanner subsystem 
24 assigns color value(s) for each of a plurality of two- 60 
dimensional sub-areas or elements of the character reference 
swatch, thereby generating an array 40 which is a color . 
mapping of the character reference material. In the illus- 
trated embodiment, scanner 24 is of a type such as that 
commercially available from Sharp Corporation as model 65 
JX-450. As explained below, array 40 comprises a plurality 
of elements 41 (corresponding to the two-dimensional sub- 
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areas), with each element 41 of array having a potential 
plurality of values assigned thereto.. 

Arrow 42 represents input of array 40 to image processor 
subsystem 26, and in particular through an input interface 50 
to image memory 52 of image processor subsystem 26 
(often abbreviated herein as "image processor"). As shown 
in FIG. 1, image processor 26 further comprises a CPU 54; 
memory 55; an input interface 58; a bus interface 59; a 
printer interface 60; a display driver 62; and, a display 
screen 64. Memory 55 includes both a ROM portion (for 
storing, e.g., executable instructions) and a RAM portion. 
An input device 70 (e.g., keyboard optionally provided with 
mouse, etc.) is connected to image processor subsystem 26 
through interface 58, Image processor subsystem 26 is 
connected through interface 60 to a printer 72. In the 
illustrated embodiment, image processor subsystem 26 is a 
textile design-type CADCAM system such as that commer- 
cially available from Sophis System. Printer 72 is preferably 
a four ink-jet type printer such an Iris model 3024 or 4012. 

FIG. 3 shows steps executed by shade range simulation 
system 20 in connection with generation of a shade range 
simulation according to a mode of the invention. Step 200 
represents a sample manufacturing run of textile material 
which, in the illustrated embodiment, is denim. While denim 
provides the example for illustration for the present inven- 
tion, it should be understood that other materials can be 
processed and simulated according to the invention, includ- 
ing non-fabrics/non-textiles. 

Step 200 of FIG. 3 typically involves production of 3,000 
to 5,000 yards of denim, or possibly as much as 20,000 yards 
of denim (in the case of a fully dye run). The denim 
production of step 200 results in a plurality of rolls of raw 
denim. At step 202, a swatch is cut from each raw roll, 
resulting in a plurality of swatches Si-S,,. At step 204, a 
standard wet process is conducted for each of the plurality 
of swatches Si-S n . 

At a random point in the production of denim occurring 
at step 200, the yam used at the random point is also utilized 
to produce two types of special test swatches. The test 
swatches of the first type, known as the WxW test swatches, 
are produced by using the random warp yarn as both the 
warp yarn and the fill yam. The test swatches of the second 
type, known as the WxF test swatches, are produced by 
using the random warp yam as the warp yam for the WxF 
test swatches and the random fill yam as the fill yam for the 
WxF test swatches. 

Generation of the WxW and WxF test swatches is 
depicted by step 206 in FIG, 3. In the illustrated embodi- 
ment, a plurality of WxW and WxF test swatches are 
generated at step 206, with each test swatch originally being 
denoted as WxW^^ or WxF ri ^ as appropriate. Since 
the test swatches are subsequentially subjected to various 
stages of processing, a test swatches is more generically 
identified as a test swatch WxW^ fl or WxF,. v (depending on 
type), wherein the subscript "s" represents a wash stage and 
the subscript "q" represents the swatch number. "Rigid" is a 
subscript which refers to a pie-wash stage of raw cut denim 
as obtained at step 206. The subscript "q" identifies a 
selected one of the plurality of test swatches. For example, 
in the illustrated embodiment, at step 206 three of each type 
of test swatches are generated: WxW WjitM , WxW^ w _ 2 , and 
WxW n .^_ 3 ; as well as WxF^, WxF_ w . 2 , and 
WxF^^.3. Each of these test swatches is subsequently 
subjected to one or more treatment stages as described 
below. 

At step 206, swatches Sj-S„ (collectively depicted as "S" 
in FIG. 2), after being washed at step 204, are each sepa- 
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rately analyzed in analyzer 20. Analyzer 20 determines 
spectral data for each swatch S h i=l, . . . n, which is then 
entered into calculation processing subsystem 22 as repre- 
sented by arrow 30 in FIG. 1. When spectral data for all 
swatches have been entered into calculation subsystem 22, 5 
at step 210 subsystem 22 averages the spectral data and, 
based on the spectral data, determines which of the swatches 
S, mostly closely corresponds to the averaged spectral data. 
At step 212, a plurality (e.g., three) of "center" or "color 
reference" swatches S R . h ^ of rigid denim are then cut 
from the same raw roll of denim that produced the swatch 
which mostly closely corresponded to the averaged spectral 
data obtained in step 210. Using this notation, the subscript 
"R" denotes a color reference swatch, the subscript "q" 
again is used to number the reference swatches. In the 
illustrated embodiment, q ranges from 1 to 3. 15 

The color reference swatches Sj^^ are subsequently 
processed in laundry or material treatment processor 28, 
along with the test swatches WxW rffid . 9 and the test 
swatches WxF^^, through one or more color-affecting 
treatment stages in the manner hereinafter described. The 
reference swatches Sx. rigb i. g and all test swatches 
WxW,.^^ and WxF,.^ are repeatedly subjected to a 
pair of sequential steps 214 and 216 (as illustrated in FIG. 
2). The number of repetitions of applications of the treat- 
ment depends on the desired treatment stage for the finished 25 
product. 

Step 214 involves a color-affecting process which is 
followed by a spectral analysis at step 216. The color- 
affecting process for the illustrated example (denim) is wet ^ 
processing (e.g., washing) of the swatches. For example, 
upon first execution, step 214 constitutes a rinse stage for 
each of the swatches S R7 WxW, and WxF. For convenience, 
after the rinse stage the swatches are referred to as S^^^ 
WxWrtnse-g* WxF rfw C - ff - After the rinse of step 214, at step 35 
216 the swatches S*.*^, WxW Wrtje . 9 , WxF rMje _ ? are each 
separately input into analyzer 20. Analyzer 20 generates 
spectral data (i.e., "primaries") for each swatch S^.^^, 
Wx ^ WxF,^,^, which spectral data is inputted to 
the calculation processing subsystem 22. 

The first execution of steps 214 and 216 can be followed 
by one or more paired executions of steps 214 and 216 for 
a second stage wash (stonewash), followed by a third stage 
wash (instant old wash), a fourth stage wash (new age wash), 
and a fifth stage wash (bleach out wash). The number of 45 
repetitions of paired steps 214 and 216 is dependent on the 
desired wash stage of the resultant product. Thus, spectral 
data ("primaries") are generated for each swatch as indicated 
by the Table 1 (and illustrated in FIG. 2). Table 1 and FIG. 
2 also show the names acquired for each of the swatches 50 
upon completion of respective wash stages. 

TABLE 1 



wash stage (s) reference watch WxW test 



WxF test 
watch mm] 



55 



rigid S R 

rinse S R 

stonewash S R 

instant old wash Sr.!^, oli ^ 
new age wash ^ 

bleach out wash Sr^ 



WxWrf^i^ 
WxWrf, 
WxW^ 
WxWfc 
WxW a< 
WxW, 



WxF* 
WxF* 
WxF„ 
WxF t 



WxR 



At step 218, the calculation processing subsystem 22 
mathematically determines, for each selected wash stage, the 
averaged color of the warp yam for the reference color 65 
swatches S^.^ of the selected wash stage. That is, for each 
wash stage, subsystem 22 determines a point in color space 



for the warp yarn for an average of the reference color 
watches S^. 7 at the end of the particular wash stage. 

As a result of step 218, for each wash stage the color for 
the warp yarn is determined at a central color candidate 
representation point in L*,a*,b* color space. See, for 
example, the color space illustration of FIG. 6. As used 
herein, the subscript "s" for point P,^ refers to washing 
stage (in the manner aforedescribed) and the subscript "0" 
distinguishes the centrally determined point from six other 
representation points (also known as elliptical or shade 
variation representation points) described below. 

In the above regard, in connection with step 218, it should 
be kept in mind that each reference swatch S*.^ is com- 
prised of both warp yam and fill yarn, with the non-colored 
fill yam having an effect on the appearance of the generally 
dyed warp yam. While it is perhaps possible that the 
reference color swatches could be dissected to extract the 
warp yarn in a manner suitable for spectral analysis, the 
mathematical processing of the invention is more advanta- 
geous. 

The determination of step 218 is accomplished by per- 
forating a spectral correction using data gleaned from the 
afore-described steps. The graphs of FIG. 5A and FIG. SB, 
together with the table of FIG. 7, represent the data utilized 
for the determination (at step 218) of the location of point 

The graphs of FIG. 5A and FIG. 5B show spectral data for 
one type of three swatches for each of the washing stages. 
FIG. 5A and FIG. SB collectively happen to illustrate 
superimposed spectral data curves for the various wash 
stages for test swatches WxF,^ (the rigid and rinse stage 
curves being shown in FIG. 5A; the instant old, new age, and 
bleachout stages being shown in FIG. SB). Although not 
necessarily visibly distinguishable for all stages, for each 
stage three curves are provided (one for each "q", q=l,2,3). 
Although illustrated only with respect to test swatches 
WxF,.^, it should be understood that similar spectral data 
curves and associated spectral data is also provided for 
reference swatches S*.^ and test swatches WxW^. 

FIG. 7 is a table showing numerical values of spectral data 
for one type of swatches (in particular for test swatches 
WxF^). Values are shown for each swatch (o=l ,23) and for 
each stage. It should be understood that a comparable table 
of data is provided for reference swatches S R . S ^ and test 
swatches WxW,. 9 . Spectral data included in such tables are 
calculated as necessary by calculating processing subsystem 
22 and stored in memory 36. 

The spectral data thus provides the calculating processing 
subsystem 22 with three sets of data— the set of data "TEST 
DATA^" illustrated in FIG. 7 for test swatches WxF a 
corresponding unillustrated set of data "REF SWATCH 
DATA^vf" for reference swatches S R t and a correspond- 
ing unillustrated set of data 'TEST DATA*^" for test 
swatches WxW . In connection with the foregoing repre- 
sentation of REF SWATCH DATA^), it should be recalled 
that the reference swatches (unlike the special test swatches 
WxW) are prepared using both warp and fill yarn, making 
the subscript "WxF* appropriate. 

Calculating processing subsystem 22 then averages the 
three values (q=l,2,3) for each stage for each type of swatch. 
So averaged, the three above-described data sets are trans- 
formed as follows: 

TEST DATAw^AVO TEST DATAuvp 
TEST DATAn^tr+AVG TEST DATA^*, 
REF SWATCH DATA^^AVG REF SWATCH DATA^ 
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Using these above-transformed sets of average data, a fur- 
ther set of averaged data (AVG REF SWATCH DATA^) 
can be mathematically obtained using conventional spread 
sheet and/or curve fitting programs, Each element of data set 
AVG REF SWATCH DATA ^ w corresponds to spectral data 
for the warp yarn required to produce the desired average 
reference swatch at the corresponding stage of washing. 
For example, a first element of set AVG REF SWATCH 
DATA^^ includes spectral data for generating the point 
P n>irf-o (i-e., swatch S/i. rtsW) ); a second element of set AVG 
REF SWATCH DATA***, (i.e., swatch 5*.^); for gen- 
erating the point P^^ol ^ so rortn - 

Thus, the determination of the color for the desired warp 
yam for each stage "s" is precisely determined as point P,_ 0 
in L*,a* t b* color space (step 218). Thereafter, a three- 
dimensional (elliptical) volume of acceptable color variation 
is determined at step 220 for a selected one of the points 
(e.g., either point P^,^. P«r ManU Mi et c). In particular, at 
step 220, calculating processing subsystem 22 obtains (e.g., 
from memory 3$) a value indicative of what degree of color 
variation from the point P,. 0 is acceptable. This acceptable 
color variation value can be prestored (e.g., in RAM portion 
of memory 55) or prompted and entered (for example) using 
input device 32. 

The elliptical volume at step 220 is essentially a deter- 
mination of DEcmc (using the spectral data for generating 
point P S _ Q ) according to conventional practice. In this regard, 
see "CMC: Calculation of Small Color Differences For 
Acceptability", AATCC Technical Manual/1992, pp. 
322-324. 

At step 222, calculation subsystem 22 determines six 
additional points in L*,a*,b* color space. These six addi- 
tional points, together with the point P^, constitute a set of 
points hereinafter referenced as color candidate representa- 
tion points. The points P^ are referred to as central color 
candidate representation points. A subset of color candidate 
representation points comprising the six additional points 
are collectively referred to as shade variation representation 
points. In particular, these additional shade variation repre- 
sentation points are chosen (as illustrated in FIG. 6) at points 
on a surface of the elliptical volume (determined at step 220) 
intersects each of secondary axes L* , ,a*' l b*' of the color 
space. The secondary axes L*'^**^*' are centered about 
point Pj. 0 and parallel to corresponding axes L*,a*,b* as 
shown in FIG. 6. The points at which the elliptical volume 
intersects the a*' axis of the color space are represented as 
P^ and ? s , 2 ; the points at which the elliptical volume 45 
intersects the b*' axis of the color space are represented as 
P,_ 3 and P^; the points at which the elliptical volume 
intersects the L*' axis of the color space are represented as 
P,. s and P,. 6 . Given the foregoing explanation (wherein 
subscript "s" is understood to refer to the selected wash 
stage), a representation point is genetically referenced as 
P,^, wherein p=0, 1 , 2, ... 6. 

At step 224, a character reference pattern is obtained. In 
the illustrated embodiment, the character reference pattern 
can be obtained in two ways. As a first way of executing step 
224, a swatch possessing the selected character of the target 
denim is scanned by scanner subsystem 24. The swatch 
scanned at step 224, identified herein as the character 
reference swatch S CR (as opposed to the color reference 
swatches S R described above), has preferably been treated 
(e.g., washed) to the same stage as the target material 
product being simulated. In an automated systems such as 
illustrated in FIG. 1, the scanning of step 224 is conducted 
in response to a scan command (carried via line 39) from the 
calculation processing subsystem 24. 

As a second way of executing step 224, rather than 
requiring a physical scan of a character reference swatch, the 
operator can select (e.g., using input device 70) which of a 
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pre-stored plurality of arrays 40 (each corresponding to 
character reference pattern) is desired for the target material- 
to-be-simulated. Thus, if desired, scanning sub-system 24 
can be eliminated or employed only selectively if the 
operator builds a library 52L of character reference swatches 
in image memory 52. 

As employed herein, the term "character" refers to one or 
more of non-color properties of the element(s), or structural 
relationship of elements, of which the target material is 
comprised. For denim, for example, character refers to one 
or more of the following: size of yam, tightness of weave, 
pattern [e.g., twill, right hand or left hand]). 

In an embodiment using scanner 24, the scanner generates 
a two-dimensional color image or color map of the scanned 
character reference swatch S CR . That is, color values are 
assigned to each of a plurality of two-dimensional sub-areas 
of the character reference swatch C CRy thereby generating a 
color array 40 for swatch S CJt , Thus, as a result of the scan , 
at step 224, each two-dimensional sub-area of the scanned 
image array is assigned shade values (i.e., an R,G,B values), 
of which there could be as many as one thousand shade 
values. Each two-dimensional sub-area is referenced herein 
as an element and preferably ultimately corresponds to a 
pixel of display screen 64. Color array 40 is stored in image . 
memory 52 of image processing subsystem 26. 

At step 228, CPU 54 of the image processing subsystem 
26 determines at what stage of washing the target material 
is to be simulated. Accordingly, at least by the time of 
execution of step 228, the image processing subsystem 26 
must be apprised as to which particular stage of washing 
(e.g., stonewash, ancient old, etc.) is to be simulated. Such 
indication of washing stage can be input through input 
device 70 (in response, for example, to a prompt and stored 
in the RAM portion of memory 55). 

At step 230 image processing subsystem 26 generates a 
plurality of swatch images SS,. p (e.g., simulated swatches) 
which represent acceptable boundary shades for the target 
denim product. As magrammatically illustrated in FIG. 6, 
each simulated swatch SS,. p corresponds to one of the color 
candidate representation points P,. p (determined at step 
222), but is corrected in view of the physical attributes of the 
character reference swatch S CR . 

Sub-steps involved in the execution of swatch simulation 
step 230 are illustrated in FIG. 8. FIGS. 10A-10E illustrate 
processing of array 40 in the course of performing step 230. 
FIG. 10A shows the status of a portion of array 40 imme- 
diately following generation of array 40 at step 226, in 
particular array elements 41 x a and 41 12 and the R,G,B 
values assigned to each element 41 1 r and 41 lr2 . 

At sub-step 230-1 (see FIG. 8), the R,G3 shade values 
obtained in step 226 are reclassified as one of a predeter- 
mined number of "reduced shades" by conventional color 
reduction techniques. In the illustrated embodiment, by way 
of non-limiting example, nine reduced shades are utilized, so 
that each display element of the scanned image array is 
associated with one of the nine reduced shades, Thus, as a 
result of the shade reduction at step 230-1, elements 41j l 
and 41 ! 2 of array 40 have reduced R.G.B values (indicated 
by R'.G^B") as seen in FIG. 10B. 

At sub-step 230-2, these nine reduced shade values are 
converted into color space (i.e., L*,a*,b* values) using 
conventional conversion techniques. Of these converted 
values, only the lightness (L*) value is utilized. Thus, 
remaining after sub-step 230-2 are nine L* reduced (sub- 
script "R") values 

!*-«(*>*• L-3(«*« • • • Lc(jr>*« Wi<h)*' • ■ • L*4<*)* 

which constitute a range R L . of lightness values centered 
about median value L^* and separated by intervals Iv as 
shown in FIG. 9, the subscript V values running from —4 to 
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44. These nine L* reduced values are generically repre- 
sented as L* w (followed by appropriate element subscript) 
in FIG. IOC. 

At sub-step 230-3, image processing subsystem 26 deter- 
mines which swatch image SS^ is first to be generated (i.e., 5 
which representation point P,. p is first to be processed). In 
the illustrated embodiment, representation point P,. 0 is first 
selected, so that swatch image SS S . Q will be generated first 
In connection with the selection of points in connection with 
sub-step 230-3 and the ensuing sub-steps, it should be 10 
recalled that the stage "s" was previously selected (see the 
description of step 228). 

For the selected representation point P^, at sub-step 
230-4, each element of array 40 in image memory 52 is 
assigned three transformed coordinate values, herein gener- 15 
ally represented as Ly*, a^*, and b* T values (see FIG. 10D). 
As explained below, the transformation of step 230-4 can 
occur either in a linear mode or in a non-linear mode, 
depending on the desired degree of precision and (in some 
applications) the available information. 20 

In both the linear mode and non-linear mode versions of 
step 230-4, only the L* values for the nine reduced shades 
are initially employed. In other words, the a* and b* values 
shown in FIG. 10C are discarded. 

In the linear mode, a range R^ of the nine lightness 25 
values is centrally re-justified about the L* coordinate value 
(i.e., Lp^p*) for a selected one of the shade variation 
representation points P,_ p (recall the elliptical representa- 
tions points were determined at step 222), That is, as 
illustrated in FIG. 9, by maintaining the intervals the 
lightness coordinates for the eight non-central reduction 
shades are re-calculated with respect to L^ p . As a result of 
the re-justification depicted by block 230-1-3, nine P-justi- 
fied or P-scaled values 
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are obtained. These nine P-justified values are generically 
referenced as L n(P ,. p , wherein n=l, . . . 9. 

Thus, in connection with the linear mode of transforma- 
tion represented by FIG. 10D, at sub-step 230-4, array 40 p 
has three values assigned to each of its elements. The first 
value (represented by subscripted values L/ in FIG. 10D) 
assigned to an element of array 40 p is an appropriate one of 
the nine P-justified or P-scaled lightness values listed above. 
The other two values (represented by subscripted values of 
a/ 1 , b T * in FIG. 10D) assigned to all elements of array 40 p 
are the a* and b* values of the selected shade variation 
representation point P,-. These non-lightness coordinate 
values of the selected shade variation representation point 
P,. p are represented as a*^ and b* Ps . Q . Thus, the a* and b* 
values are constant for each element of array 40 p . That is, for 
all elements in array 40 p , z T *=a* Ps _ p and b T *=b* Fj . p . The 
lightness correction implemented by this assignment is thus 
dependent upon the lightness of the character reference 
swatch S Cft . 

In the non-linear mode of executing sub-step 230-4, nine 
L* reduced (subscript "R") values 

from step 230-2 are utilized. As an example of the non-linear 
mode of execution of sub-step 230-4, consider the central 
color candidate representation point P,^ illustrated in FIG. 
6. The L*, a*, and b* coordinate values of point P,. 0 are d5 
known as a result of execution of step 218 of FIG. 3. In the 
non-linear mode, those elements of array 40 having a 



lightness value which is the central lightness value L^* 
acquire the same a* and b* coordinate values as central color 
candidate representation point P,.^ The elements 41 of array 
40 having lightness values other than Lq {J{ * are assigned 
other a*, b* values as described below. 

In the above regard, in addition to the calculation results 
obtained from step 218 (see FIG. 3), the image processor 26 
has access to the spectral data collected by calculation 
processing subsystem 24 at step 216. This spectral data, 
having exemplary values such as shown in FIG. 7, forms a 
plurality of curves such as illustrated in FIG. 5A and FIG. 
5B. One of the curves corresponds to the point P,_ 0 calcu- 
lated at step 218. In connection with the non-linear mode, 
the curve corresponding to point P,_ 0 and eight other points/ 
curves are utilized. In particular, taking the curve corre- 
sponding to point ? s . 0 as central point, the four most 
neighborly curves on each side of the curve corresponding 
to point P^o are taken as neighboring curves (e.g., four 
curves below and four curves above). 

FIG. 11 illustrates further sub-steps involved in the non- 
linear execution of sub-step 230-4. In particular, at sub-step 
230-4-1-NL, the eight neighboring curves are determined as 
described above. Then, at sub-step 230-4-2-NL, a conven- 
tional K/S calculation is made with respect to each of the 
eight neighboring curves. 

At sub-step 230-4-3-NL, image processor 26 determines 
which of the color candidate representation points is next to 
be processed. If it is determined at sub-step 230-4-4-NL that 
the next point is the central point P s . 0 , then sub-steps 
230-4-5-NL through 230-4-7-NL are executed. Otherwise, 
sub-steps 230-4-8-NL through sub-step 230-4-11-NL are 
executed. 

At sub-step 230-4-5-NL, image processor 26 determines 
which of the eight non-central L* reduced (subscript "R") 
values needs to have a*, b* values assigned thereto. For 
example, for a first execution of sub-step 230-4-5-NL, 
elements having the lightness value L_ ltf?) * may be pro- 
cessed. 

The actual assignment of non-L* coordinate values per 
element occurs at sub-step 230-4-6-NL. At sub-step 230-4- 
6-NL, image processor 26 determines from the set of nine 
K/S-converted curves, and the associated L*, a*, and b* 
values, what a*, b* values should be assigned to array 
elements having the lightness currently being processed. For 
45 example, the a* and b* values of the central point Pj. 0 have 
been previously calculated (see step 218 in FIG. 3). If the 
curves do not having exactly corresponding L* values for 
the L„ (/ > x _ p) being processed, an interpolation can be per- 
formed according to standard interpolation techniques. 

In connection with the foregoing, it is advantageous to 
conduct sub-step 230-4-6-NL in terms of an L, c and h 
(lightness, chroma, and hue) polar coordinate system rather 
than the L*, a*, b* coordinate system. In this regard, it is 
well -recognized thai L*, a*, and b* values can also be 
represented in another (polar) coordinate system which has 
the axes L*, c (chroma) and h (hue). Accordingly, it will be 
understood that the man skilled in the art can make the 
necessary conversions to L, c, h coordinates, so that (for the 
most part) the remainder of the non-linear mode will be 
described with respect to the L, c, h coordinate system 

After an execution of sub-step 230-4-6-NL, image pro- 
cessor 26 checks at sub-step 230-4-7-NL to determine 
whether all non-central L* reduced (subscript "R") values 
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55 



60 



L-^^)*, L. 



•-3<*)*> ■ • • L-k*)*. . . . L^^* have been 

processed. If not, execution loops back to sub-step 230-4- 
5-NL for acquisition of yet another non-central L* reduced 
(subscript "R") value. If all elements have been processed, 
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execution returns to sub-step 230-4-3-NL for obtaining yet 
another color candidate representation point. 

As indicated above, when a non-central color candidate 
representation point is processed, sub-steps 230-4-8-NL 
through sub-step 230-4- 12-NL are executed for a shade 5 
variation representation point. Such shade representation 
points are represented as points P,.p, p=l, . . . 6 in FIG. 6. 

At sub-step 230-4-8-NL, image processor 26 obtains the 
c and h values for the previously calculated (see step 222 of 
FIG. 3) shade variation representation point P,_p currently 10 
being processed Further, at sub-step 230-4-8-NL, image 
processor 26 assigns all elements 41 of array 40 having the 
central L* reduced (subscript "R") value the c and h values 
of the shade variation representation point P,_ p currently 
being processed. 15 

At sub-step 230-4-9-NL, the spectral curves previously 
employed (e.g., at sub-step 230-4-6-NL) are adjusted using 
a coordinate offset factor. The coordinate offset factor uti- 
lizes a positional differential between the particular shade 
variation representation point being processed P^ p and the 20 
central color candidate representation point P,^. 

In the above regard, consider a shade variation represen- 
tation point P,.^ and a central color candidate representation 
point P^ 0 as in FIG. 13. FIG. 13 shows an ellipse, such as 
the ellipse of FIG. 6, but in only two dimensions (c, h) in a 25 
L-c-h coordinate system. As illustrated in FIG. 13, shade 
variation representation point P^ and central color candi- 
date representation point P,^ have essentially the same 
lightness and hue, so that the coordinate offset factor is the 
Ac. Accordingly, for the "c" data used to generate the eight 30 
neighboring curves with respect to shade variation repre- 
sentation point ? s , x is adjusted or offset by the "c" coordinate 
offset factor. For shade variation representation point P^, 
Ah and AL are zero. 

In connection with the adjustment of curves at 230-4-9- 35 
NL, it should be understood that a similar coordinate offset 
factor can be determined and used for other shade variation 
representation points. For example, for shade variation rep- 
resentation point P,.y, a coordinate offset factor Ah would be 
determined and used to adjust the "h" values used for the 40 
spectra curves. 

At sub-step 230-4-10-L, much like sub-step 230-4-5-NL, 
image processor 26 determines which of the nine L* reduced 
(subscript "R") values of array 40 need next be processed. 
For example, for a first execution of 230-4-10-NL, perhaps 45 
the non-central L* reduced (subscript "R") value L_ 1(/J) is 
processed. For each non-central L n(/0 *, sub-steps 230-4-11- 
NL is executed. 

At sub-step 230-4- 11 -NL, image processor 26 makes an 
actual assignment of c and h values to the elements 41 of so 
array having the L n(J?) value being processed (for example, 
the non-central L* reduced (subscript "R") value L_ 1(/r) is 
processed for the first execution of sub-step 230-4-11-NL). 
The assignments of sub-step 230-4- U-NL are similar to 
those of sub-step 230-4-6-NL, it being understood however 55 
that the adjusted spectral curves (adjusted as in sub-step 
230-4-9-NL) are utilized. 

At sub-step 230-4-12-NL, image processor 26 checks to 
determine if all non-central L* reduced (subscript "R") 
value L n(R) have been processed. If not, execution returns to 60 
sub-step 230-4-10-NL to obtain yet another non-central L* 
reduced value. If all lightness values have been processed 
for the current shade variation representation point, process- 
ing returns to sub-step sub-step 230-4-3-NL to obtain 
another shade variation representation point 65 

Thus, in the non-linear mode of executing sub-step 230-4 
as described above, the L w * values for the nine reduced 
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shades (see FIG. 10C and sub-step 230-2) are utilized as the 
Lr* values (see FIG. 10D). However, the a r * and b r * values 
are obtained from spectral data previously generated at step 
216 (see FIG. 3), and which is interpolated and adjusted as 
necessary. 

In connection with the non-linear mode of execution of 
sub-step 230-4, particularly sub-steps 230-4-6-NL and 230- 
4-11-NL, ordinarily neighboring curves are utilized. For 
example, if the stage to be simulated is stonewash, it is likely 
that coordinate values not only from the stonewash spectral 
data, but also from the rinse and instant old spectral data will 
be utilized. For extreme wash stages (e.g., bleachout), coor- 
dinate values from the new age stage and other extrapolated 
values will be utilized. 

Having just described both the linear and non-linear 
modes of executing sub-step 230-4, discussion returns to the 
remaining sub-steps of FIG. 3. 

At subrstep 230-5, each display element of desired image 
simulation array 40 (in color space values) is converted back 
into R,G,B values, yielding array 40p^ Gfl , Such conversion 
is basically the inverse of the conversion of sub-step 230-2. 
Conversion values R^G^ are illustrated in FIG. 10E. 

At sub-step 230-6, the array 40p.^ Cfl is visibly reproduced 
as simulated swatch SS^ p using either or both of printer 72 
or display screen 64. Array 40^^ is replicated on display 
screen 64 (using display driver 62) and on (e.g., paper) 
hardcopy (by printer 72) in conventional manner. 

Upon replication of a simulated swatch, at step 230-8 the 
image processing subsystem 26 determines if simulation 
swatches have been generated for all seven representation 
points. If simulation swatches have yet to be generated for 
other representation points, the next representation point is 
selected at step 230-8. For example, after representation 
point P^o has been processed, representation point P,., 
might next be processed, followed by point P,_ 2 , and so 
forth, 

Hence, sub-steps 230-4 through 230-8 constitute a loop, 
most of which is executed for each representation point in 
connection with the generation of an associated simulated 
swatch SSj_ p , p=0, 1, ... 6. After all representation points 
have been thusly processed as determined at sub-step 230-7, 
generation step 230 is terminated (as indicated by sub-step 
230-9). 

Upon completion of step 230, image processing sub- 
system 26 determines if simulated swatches are needed for 
another stage (step 232). If so, execution returns to step 228 
for a determination of what stage should next be processed. 
Thus, it should be understood from the foregoing that a set 
(e.g., seven) of simulated swatches can be generated for one 
or more, .or even all, washing stages. When all swatch 
generation is complete, the process of FIG. 3 terminates 
(step 234). 

It should also be understood that the procedure described 
above can be further automated to provide, for example, for 
the subsequent scanning of other character reference 
swatches. 

Thus, upon completion of step 230, the operator is pro- 
vided — prior to manufacture — with a set of simulated 
swatches SS_ p which visibly depict acceptable boundary 
shades for the target denim product 

Using the central simulated swatch SS^ and the bound- 
ary shade simulated swatches SS^p, p=l, ... 6, a customer 
can make an accurate deterrnination of shade acceptability 
even before rnanufacture. Similarly, quality control person 
can have a more comprehensive frame of reference in which 
to determine if production material satisfies the expectation 
of the customer. 
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Whereas the foregoing discussion concerns the generation 
of a colored target material using judiciously selected color 
candidate representation points, it should be understood that 
the invention can also generate simulations based on other 
criteria. For example, as shown in FIG. 3B, in an alternate 5 
embodiment, formula performance data gleaned from a 
color formula analysis system (such as that marketed by 
SheLyn, Incorporated) can instead be input at step 322 to 
image processing sub-system 26 as formula performance 
representation points. The formula performance data 
includes data represented by the display of FIG. 11. That is, 
in order to generate a set of swatches at step 330, a plurality 
of representation points in color space are input at step 322 
based on the data set that produced the output shown in FIG. 
U. For example, the representation points input at step 322 
might include the lightest and darkest points from the output 15 
data set of FIG. 11 as well as points having the extreme-most 
readings along the a* and b* axes (in both the negative and 
positive sense). Thus, a total of six representation points 
would be selected, such that a set of six formula perfor- 
mance simulated swatches would be produced at step 330. 20 

In connection with the steps executed in FIG. 3, it should 
be understood that execution of steps 324, 326, 328, and 334 
is understood from the previous description of the execution 
of steps 224, 226, 228, and 234, respectively. Sub-steps 
encompassed within step 330 are similar to those illustrated 25 
(in FIG. 8) for step 230. 

While the invention has been particularly shown and 
described with reference to the preferred embodiments 
thereof, it will be understood by those skilled in the art that 
various alterations in form and detail may be made therein 3° 
without departing from the spirit and scope of the invention. 

The embodiments of the invention in which an exclusive 
property or privilege is claimed are defined as follows: 

1. A method for simulating a colored target material, the 
method comprising: 35 

selecting a point in color space, the selected point having 

at least one coordinate value; 
storing in a memory character reference values for a 

character reference pattern having at least one desired ^ 

character property of the target material; 
using the coordinate value of the selected point to modify 

the reference values stored in the memory; 
using the reference values stored in the memory to 

provide a visual output representative of the target 45 

material. 

2. The method of claim 1, wherein the selected point has 
a lightness coordinate L* in CUE color space and non- 
lightness coordinate values a*, b* in CBB color space. 

3. The method of claim 1, wherein using the reference 50 
values stored in the memory to provide a visual output 
comprises printing the reference values stored in the 
memory to provide a simulation of the colored target mate- 
rial. 

4. The method of claim 1, wherein using the reference 55 
values stored in the memory to provide a visual output 
comprises driving a screen display to provide a simulation of 
the colored target material. 

5. The method of claim 1, wherein the step of obtaining 
the character reference pattern comprises analyzing a char- 60 
acter reference material by mapping the character reference 
material into an array and, based on character properties of 
the character reference material, assigning at least initial 
lightness values to each of a plurality of elements of the 
array. 65 

6. The method of claim 1, wherein the step of obtaining 
the character reference pattern comprises selecting a desired 



character reference pattern from a plurality of character 
reference patterns stored in a memory. 

7. The method of claim 1, wherein the point in space is 
selected from a set of points representing potential shade 
variations of a dye formula for the target material. 

8. The method of claim 7, wherein each point in the set of 
points has an L* coordinate value, an a* coordinate value, 
and a b' coordinate value, and wherein the selected point is 
one of the following points: 

the point in the set with the lightest coordinate value of 
potential shade variation of the dye formula; 

the point in the set with the darkest coordinate value of 
potential shade variation of the dye formula; 

the point in the set with the greatest a* coordinate value; 

the point in the set with the least a* coordinate value; 

the point in the set with the greatest b* coordinate value; 

the point in the set with the least b* coordinate value. 

9. The method of claim 1, wherein the colored target 
material is a desired stage of finish of denim. 

10. A method for simulating a colored material, the 
method comprising: 

storing, in a memory, character reference values repre- 
senting a character reference pattern having at least one 
desired character property of the material-to-be-simu- 
lated; 

using a coordinate value of a selected point in color space 
to modify selected ones of the character reference 
values in the memory and thereby yield values for a 
modified pattern; 

using the values of the modified pattern to provide a visual 
output representative of the material-to-be-simulated. 

11. The method of claim 10, wherein the selected point in 
color space is obtained by conducting spectral analysis for 
material subjected to at least one stage of a color-affecting 
process. 

12. The method of claim 11, wherein the material is fabric, 
and wherein the selected point in color space is obtained by 
conducting spectral analysis for fabric subjected to at least 
one stage of a washing process. 

13. The method of claim 12, wherein the fabric is denim. 

14. A method for simulating a range of shade variations 
for a colored target material, the method comprising: 

selecting a point in color space corresponding to a desired 
constituent color of the colored target material; 

determining an acceptable color variation about the target 
point; 

selecting a plurality of shade variation representation 
points in color space in accordance with the acceptable 
color variation, thereby forming a set of color candidate 
representation points comprising the selected point and 
the shade variation representation points; 

generating a simulation output for each representation 
point, thereby forming a set of simulation outputs 
depicting a range of predetermined shade variation. 

15. The method of claim 14, wherein the color space has 
L*, a*, and b* axes and wherein the shade variation repre- 
sentation points lie at intersections of the three dimensional 
volume with a secondary axis of color space, the secondary 
axis being parallel to one of the L*, a*, b* and axes in CLE 
color space. 

16. The method of claim 14, wherein the point has both 
a lightness coordinate value and non-h'ghtness coordinate 
values, further including 

obtaining a character reference pattern having at least one 
desired character property of the target material, the 
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character reference pattern being stored as an array, 
each of a plurality of elements of the array being 
assigned a corresponding initial lightness value; 
and wherein, to form the set of simulation outputs depict- 
ing a range of predetermined shade variation, the 5 
following steps are performed for each color candidate 
representation point: 

(a) using the lightness coordinate value of the color 
candidate representation point to form adjusted light- 
ness values for each of the elements of the array; 10 

(b) transforming the array by assigning the non-light- 
ness coordinate values of the color candidate repre- 
sentation point to each element of the array whereby, 
together with the adjusted lightness values the non- 
lightness coordinate values of the representation 
point, each element of the transformed array is 1 
representable in color space; and 

(c) using the transformed array to provide a visual 
output representative of a corresponding color shade 
variation of the target material. 

17. The method of claim 14, wherein the target point has 20 
both a lightness coordinate value and non-lightness coordi- 
nate values, the selected point belonging to a predefined set 

of points, the method further including: 

obtaining a character reference pattern having at least one 
desired character property of the target material, the 25 
character reference pattern being stored as an array, 
each of a plurality of elements of the array being 
assigned a corresponding initial lightness value; 

and wherein, to form the set of simulation outputs depict- 
ing a range of predetermined shade variation, the 30 
following steps are performed for each representation 
point: 

(a) for a selected element of the array, using the lightness 
coordinate value of the selected element to determine 
transformation non-lightness coordinate values to be 
assigned to the selected element, the transformation 
non-lightness coordinate values to be assigned to the 
selected element being determined by referencing the 
predefined set of points in color space, the determina- w 
tion utilizing, as the transformation non-lightness coor- 
dinate values for the selected element, non-lightness 
coordinate values of a matched one the points in the set, 
the matched point in the set having a lightness coordi- 
nate value which fulfills a predetermined criteria with 
respect to the lightness coordinate value of the selected 
element; 

(b) transfonning the array by assigning the non-lightness 
coordinate values of the representation point to the 
selected element of the array whereby, together with the 50 
lightness cxwrdinate values of the selected element of 
the array, the selected element of the transformed array 

is representable in color space; and 

(c) using the array with the transformation non-lightness 
coordinate values for the selected element to provide a 55 
visual output representative of the target material. 

18. The method of claim 17, wherein the predetermined 
criteria is closest numerical proximity to the lightness coor- 
dinate value of the selected element 

19. The method of claim 17, wherein the predefined set of 60 
points in color space is obtained by conducting spectral 
analysis for material subjected to at least one stage of a 
color-affecting process. 

20. The method of claim 19, wherein the material is 
fabric, and wherein the predefined set of points in color 65 
space is obtained by conducting spectral analysis for fabric 
subjected to at least one stage of a washing process. 
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21. The method of claim 20, wherein the fabric is denim. 

22. A method for simulating a colored target material 
expected to result from at least one color-affecting treatment, 
the target material being composed of a first constituent 
element and a second constituent element, the method 
comprising: 

(a) analyzing a colored reference material, the colored 
reference material being composed of the first constitu- 
ent element and the second constituent element; 

(b) analyzing a first test material, the first test material 
being composed only of the first constituent element; 

(c) analyzing a second test material, the second test 
material being composed of the first constituent ele- 
ment and the second constituent element, the first 
constituent element of the second test material being 
the same as that of the first test material; 

(d) subjecting the color reference material, the first test 
material, and the second test material to at least one 
color-affecting treatment; 

(e) analyzing each of the color reference material, the first 
test material, and the second test material after the at 
leash one color-affecting treatment; 

(f) using analyses from steps (a)-(c) and (e) to determine 
a color which corresponds to the color for the first, 
constituent element, after the at least one color-affect- 
ing treatment. 

23. The method of claim 22, wherein the target material 
is a textile, and wherein the first constituent element is warp 
yam and the second constituent element is fill yarn. 

24. The method of claim 22, wherein a plurality of color 
reference materials, a plurality of first test materials, and a 
plurality of second test materials are utilized. 

25. The method of claim 22, wherein a plurality of 
color-affecting treatments are preformed, and wherein step 
(e) is performed after each color-affecting treatment. 

26. The method of claim 25, wherein the target material 
is denim, wherein the first constituent element is warp yarn 
and the second constituent element is fill yam, and wherein 
the color-affecting treatments are sequential washing pro- 
cesses. 

27. The method of claim 25, wherein the selected point 
has both a lightness coordinate value and non-lightness 
coordinate values, and wherein the target material is simu- 
lated by executing the following steps: 

obtaining a character reference pattern having at least one 
desired character property of the target material, the 
character reference pattern being stored as an array, 
each of a plurality of elements of the array being 
assigned a corresponding initial lightness value; 

using the lightness coordinate value of the selected point 
to form adjusted lightness values for each of the 
elements of the array; 

transforming the array by assigning the non-lightness 
coordinate values of the selected point to each element 
of the array whereby, together with the adjusted light- 
ness values and the non-lightness coordinate values of 
the selected point, each element of the transformed 
array is representable in color space; and 

using the transformed array to provide a visual output 
representative of the color of the target material. . 

28. Apparatus for simulating a colored target material, the 
apparatus comprising: 

a memory wherein is stored character reference values for 
a character reference pattern, the character reference 
pattern having at least one desired character property of 
the target material; 
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an image processor which uses a coordinate value of the 
a selected point in color space to adjust at least some of 
the character reference values and thereby form an 
adjusted character pattern; 

an output device which uses the adjusted character pattern 5 
to provide a visual output representative of the target 
material. 

29. The apparatus of claim 28, wherein the selected point 
has a lightness coordinate L* in CIE color space and 
non-lightness coordinate values a*, b* in CIE color space. 10 

30. The apparatus of claim 28, wherein the output device 
is a printer. 

31. The apparatus of claim 28, wherein using the output 
device is a display screen. 

32. The apparatus of claim 28, further comprising a 
scanner which analyzes the character reference material by 15 
mapping the character reference material into an array and, 
based on character properties of the character reference 
material, assigns at least initial lightness values to each of a 
plurality of elements of the array. 

33. Apparatus for simulating a colored target material, the 20 
apparatus comprising: 

a first memory wherein is stored values for a character 
reference pattern having at least one desired character 
property of the target material; ^ 

a second memory wherein is stored at least one coordinate 
value for a selected point in color space; 

an image processor which uses the selected point in color 
space to adjust at least some of the stored values for the 
character reference pattern, thereby yielding an 30 
adjusted pattern; and 

an output device which uses the adjusted pattern to 
provide a visual output representative of the target 
material. 

34. The apparatus of claim 23, wherein the output device 35 
is a printer. 

35. The apparatus of claim 33, wherein using the output 
device is a display screen, 

36. The apparatus of claim 33, further comprising a 
scanner which analyzes a character reference material by 40 
mapping the character reference material into an array in the 
memory and, based on character properties of the character 
reference material, assigns at least initial lightness values to 
each of a plurality of elements of the array. 

37. Apparatus for simulating a colored target material, the 45 
apparatus comprising: 

a memory for storing therein at least one coordinate value 
for a selected point in color space corresponding to a 
desired color of the colored target material; 

an input device for inputting criteria for developing an 50 
acceptable color variation about the selected point; 

a processor which uses the inputted criteria to further 
select a plurality of shade variation representation 
points in color space, thereby forming a set of rcpre- 55 
sentation points comprising the selected point and the 
plurality of shade variation representation points; and 

an output device for generating a simulation output for 
each representation point, thereby forming a set of 
simulation outputs depicting a range of predetermined 60 
shade variation. 

38. Hie apparatus of claim 37, wherein the output device 
is a printer. 

39. The apparatus of claim 37, wherein using the output 
device is a display screen. 65 

40. A method for simulating a colored target material 
expected to result form at least one color-affecting treatment, 
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the target material comprised of a first constituent element 
and a second constituent element, the method comprising: 

(a) analyzing a plurality of test materials, including a first 
test material and a second test material, the first test 
material comprising the first constituent element but 
not the second constituent element, the second test 
material comprising the second constituent element but 
not the first constituent element, the first constituent 
element of the second test material being the same as 
that of the first test material; 

(b) subjecting the color reference material, the first test 
material, and the second test material to at least one 
color-affecting treatment; 

(c) analyzing the first test material and the second test 
material after the at least one color-affecting treatment; 

(d) using the analyses from steps (a) and (c) to determine 
a color which corresponds to the color for the first 
constituent element after the at least one color-affecting 
treatment. 

41. The method of claim 40, wherein the target material 
is a textile, and wherein the first constituent element is warp 
yam and the second constituent element is fill yam 

42. The method of claim 40, wherein a plurality of color 
reference materials, a plurality of first test materials, and a 
plurality of second test materials are utilized. 

43. The method of claim 40, wherein the target material 
is denim, wherein the first constituent element is warp yam 
and the second constituent element is fill yam, and wherein 
the color-affecting treatments are sequential washing pro- 
cesses. 

44. A method for simulating a color target material, the 
method comprising: 

selecting a point in color space and generating signals 
therefor; 

obtaining signals for a character reference pattern, the 
character reference pattern having at least one desired 
character property of the target material; 

using a processor to generate simulation signals which 
simulate the colored target material, the processor 
generating the simulation signals by adjusting at least 
some of the signals for the character reference pattern 
in accordance with the signals for the point in color 
space; and 

driving an output device with the simulation signals to 
obtain a visually-perceptible image of the colored tar- 
get material. 

45. A method for simulating a range of shade variations 
for a colored target material, the method comprising: 

(1) inputting into a processor signals indicative of a 
selected point in color space; 

(2) inputting into the processor signals indicative of an 
acceptable color variation about the selected point; 

(3) using the processor to generate, from the signals input 
in steps (1) and (2), signals for a plurality of shade 
variation representation points in color space, thereby 
generating a set of signals for a corresponding set of 
color candidate representation points comprising the 
selected point and the shade variation representation 
points; 

(4) using the set of signals generated at step (3) to drive 
an output device and thereby provide a corresponding 
set of visually-perceptible simulation outputs for each 
representation point, and thereby depicting a range of 
predetermined shade variation. 

***** 
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